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ABSTRACT 


A  new  approach  for  solving  electromagnetic  wave  propagation  problems  is  currently 
being  developed  at  the  Naval  Command,  Control  and  Ocean  Surveillance  Center 
(NCCOSC),  RDT&E  Division  (NRaD).  This  new  approach  is  based  upon  an  8  by  8  matrix 
representation  of  the  Maxwell  field  equations.  In  addition,  a  computer  software  package 
based  on  this  matrix  representation  of  electromagnetic  theory  is  also  being  written  and 
tested  at  NRaD  to  handle  a  variety  of  scenarios  involving  electromagnetic  wave 
propagation  through  matter.  This  software  package  is  referred  to  as  the  MATURE 
Program.  MATURE  is  the  acronym  for  Matrix  Approach  To  Understanding  Relativistic 
Electrodynamics.  The  MATURE  Program  is  written  in  MATLAB  code  for  use  on  a  Sun  4 
SPARCstation  2  workstation.  Under  Independent  Research  (IR)  FY  92  funding,  this 
matrix  approach  was  successfully  employed  in  solving  problems  dealing  with 
electromagnetic  wave  propagation  through  dielectric,  crystalline,  linear  electro-optic, 
and  magneto-optic  materials  of  infinite  extent.  Under  the  Office  of  Naval  Research 
(ONR)  FY  93  funding,  this  matrix  formulation  was  extended  to  handle  problems 
involving  wave  propagation  through  multilayer  dielectric  media  with  planar 
boundaries.  Presented  in  this  technical  document  is  the  underlying  theory  of  this  matrix 
approach.  Several  numerical  examples,  based  on  the  use  of  the  MATURE  Program,  are 
also  included  to  illustrate  the  use  of  the  matrix  approach  in  solving  electromagnetic 
wave  propagation  problems. 
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1.0  INTRODUCTION 


In  Part  I  of  this  series  (Bocker  &  Frieden,  1992),  a  new  matrix  formulation  of  classical 
electromagnetic  theory  for  vacuum  is  presented.  The  basis  of  this  formulation  is  a  skew- 
Hermitian  space-time  8  by  8  differential  matrix  operator.  There  it  was  shown  that  the 
vector  form  of  the  Maxwell  field  equations  for  vacuum  can  be  cast  into  a  compact  matrix 
form.  From  the  matrix  form  of  the  Maxwell  field  equations,  other  fundamental  results  of 
electromagnetic  theory  are  easily  derived  with  use  of  the  simple  matrix  multiply 
operation:  (1)  the  electromagnetic  wave  and  charge  continuity  equations;  (2)  the  Lorentz- 
and  Coulomb-gauge  definitions  of  the  electromagnetic  potentials;  (3)  the  wave  equations 
for  the  potentials;  and  (4)  Poynting's  theorem  on  energy  conservation.  Taking  the  four¬ 
dimensional  Fourier  transform  of  the  matrix  form  of  Maxwell's  equations  leads  to:  (5)  a 
Fourier  representation  of  the  Maxwell  field  equations;  (6)  their  inversion,  for  the  fields 
directly  in  terms  of  the  sources  in  Fourier-space;  and  (7)  corresponding  inversion 
formulae  in  direct-space  through  the  use  of  the  convolution  theorem.  In  Part  II  of  this 
series  (Bocker,  1992),  the  8  by  8  electromagnetic  matrix  formulation  for  vacuum  was 
extended  to  include  the  presence  of  matter.  Emphasis  was  placed  on  electromagnetic 
wave  propagation  through  linear,  homogeneous,  and  anisotropic  optical  materials  of 
infinite  extent  without  boundaries.  The  matrix  form  of  the  Maxwell  field  equations  was 
cast  into  an  8  by  8  matrix  eigenvalue  representation.  In  this  technical  document,  we 
present  the  underlying  theory  for  solving  electromagnetic  wave  propagation  problems 
involving  multilayer  dielectric  materials  with  parallel  planar  boundaries.  The  matrix 
representation  of  the  Maxwell  field  equations  plays  a  central  role  in  this  investigation  as 
well.  As  will  be  shown,  boundary  conditions  are  automatically  taken  into  account  by 
partitioning  both  the  electromagnetic  wavefunctions  and  the  dielectric,  permeability, 
and  conductivity  matrices  of  the  dielectric  medium  using  Heaviside  unit  step  functions. 

Because  of  the  power  of  matrix  operations,  many  of  the  fundamental  relationships 
and  results  in  electromagnetic  wave  propagation  theory  are  easily  derived  without  the 
need  for  the  usual  plethora  of  vector  calculus  identities  that  have  become  the  standard  in 
these  derivations.  Instead,  mathematical  operations  involving  matrix  multiplication, 
matrix  inversion,  and  eigensolutions  of  matrices  can  be  employed.  Available  off-the- 
shelf  computer  software  packages,  like  MATLAB  and  MATHEMATICA,  are  well  suited 
for  these  matrix  operations.  A  computer  software  package  entitled  MATURE  is  currently 
being  written  and  tested  at  NRaD  to  handle  a  variety  of  scenarios  involving 
electromagnetic  wave  propagation  through  matter.  MATURE  is  the  acronvm  for  Matrix 
Approach  To  Understanding  Relativistic  Electrodynamics.  The  MATURE  program  is 
being  written  in  MATLAB  code  for  use  on  a  Sun  4  SPARCstation  2  workstation. 
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2.0  THEORY 


In  this  section,  we  will  derive  the  equations  for  predicting  the  reflection  and 
transmission  characteristics  of  a  multilayer  dielectric  medium  with  planar  boundaries 
illuminated  by  a  monochromatic  electromagnetic  plane-wave.  We  will  begin  with  the 
vector  form  of  the  Maxwell  field  equations  for  matter.  The  vector  form  of  the  Maxwell 
field  equations  will  then  be  cast  into  the  compact  matrix  form  described  in  Parts  1  and  II 
of  this  series.  The  matrix  form  will  then  be  specialized  to  dielectric  materials.  Next,  the 
dielectric  medium  will  be  partitioned  into  a  multilayer  medium  with  planar  boundaries 
using  Heaviside  unit  step  functions.  As  will  be  shown,  the  electromagnetic  boundary 
conditions  at  each  of  the  planar  interfaces  will  be  automatically  taken  into  account  using 
this  approach.  Then,  monochromatic  electromagnetic  plane-wave  solutions  will  be 
considered.  This  will  lead  to  system  operator  description  of  wave  propagation  through 
the  multilayer  medium.  The  operator  description  will  then  be  used  to  calculate  the 
reflection  and  transmission  characteristics  of  the  multilayer  medium  for  arbitrarily 

-  J 

chosen  wavelengths,  angles  of  incidence,  and  polarization  states  of  the  incident 
electromagnetic  radiation. 

2.1  THE  MAXWELL  FIELD  EQUATIONS 

The  fundamental  equations  of  classical  electromagnetic  phenomena,  namely  the 
Maxwell  field  equations,  serve  as  our  starting  point.  In  the  Gaussian  system  of  units,  the 
four  Maxwell  field  equations  in  vector  form  are  given  by  the  following  (Jackson,  1962): 

Ampere-Maxwell  law 

VxH(r.f)  =  -|rD(r,f)  +—Je(r,t)  (1) 

cat  c 

Gauss'  law  for  electricity 

V*D(r,f)  =4npe(r,t)  (2) 

Faraday's  law  of  induction 

Vx  E(r,t)  =  -  - ~~B  (r,t)  ~—fm(r,t)  (3) 

c  at  c 

Gauss'  law  for  magnetism 


V«B  (r,  t)  =  4npm  (r,  t) . 


(4) 


The  physical  quantities  appearing  in  the  Maxwell  field  equations  are:  E(r.t)  the 
electric  field;  D  (r,  t)  the  electric  displacement;  B  ( r.  t)  the  magnetic  induction;  H(r.t) 
the  magnetic  field;  Je  ( r,  t )  the  electric  current  density;  /'"  ( r,  / )  the  magnetic  current 
density;  p e{r,t)  the  electric  charge  density;  p"'(r.  t)  the  magnetic  charge  density; 
(r,  t)  a  space— time  point;  V  the  del  operator;  and  c  represents  the  speed  of  light  in 
vacuum.  Both  magnetic  charge  and  current  densities  (Magid,  1972)  have  been  included 
in  Maxwell's  equations  for  purposes  of  completeness.  They,  of  course,  can  be  set  equal  to 
zero  since  magnetic  charge  has  not  been  discovered  in  nature.  The  electric  displacement 
and  electric  field  vectors,  as  well  as  the  magnetic  induction  and  magnetic  field  vectors, 
are  related  (Jackson,  1962)  through  the  mathematical  expressions 

D(r,t)  =  E(r,t)  +4n  P(r,t)  (5) 

and 

B  (r,  t)  =  H  (r,  t)  +  4n  M  (r,  t) ,  (6) 

where  P(r,  t)  and  M  (r,  t )  are  the  macroscopic  polarization  and  magnetization  vectors. 

The  Maxwell  field  equations  (1)  through  (4),  with  the  use  of  equations  (5)  and  (6),  can 
be  cast  (Bocker,  1992)  into  the  following  covariant  matrix  equation 


M,  M2 
M2  M, 


/i 


+  4ji 

M, 

o' 

O 

d2 

4tv 

c 


(7) 


This  matrix  representation  of  the  Maxwell  field  equations  forms  the  basis  for  all  other 
work  described  in  this  report.  For  the  Lorentz  gauge  choice  (Bocker  &  Frieden,  1992),  the 
4  by  4  differential  matrix  operators  1MJ  and  fM2]  are  defined  by  the  equation 
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where 


r  -  ( x ,  y.  z)  and  r  =  ict. 


(9) 


The  imaginary  quantity  >  is  equal  to  y-l .  The  matrix  |  0)  in  equation  (7)  represents  the 
4  bv  4  null  matrix.  The  4  by  1  electric  and  magnetic  field  vectors  i  f,  |  and  |  f:l  appearing 
in  equation  (7)  are  defined  by 


r  i 

£* 

E 

and  [/,]  = 

H 

y 

V 

E, 

!i  -J 

h\ 

0_ 

_0_ 

(10) 


The  4  by  1  polarization  and  magnetization  vectors  |<f,|  and  |  d2\  are  defined  by 


X 

Mx 

Yv 

si\ 

Py 

and 

M- 

My 

p. 

Mz 

Oj 

L 

(11) 


The  electric  and  magnetic  source  vectors  fs,]  and  [s2|  appearing  in  equation  (7)  are 
relativistic  4-vectors  defined  by 


n 


(12) 


It  is  through  the  polarization,  magnetization,  and  source  vectors  that  the  material 
characteristics  of  the  medium  are  mathematically  described.  As  previously  indicated, 
our  interest  concerns  dielectric,  i.e.,  linear,  homogeneous,  isotropic,  and  conductive 
materials  only. 
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2.2  DIELECTRIC  MEDIUM 


For  a  linear  medium,  the  polari/  .v>n  and  magnetization  \.  ect(.>rs,  |  A]  |  and  i ;/:  j ,  are 
related  to  the  electric  and  magrvi  _  held  vectors,  [/",  |  and  1 | ,  by  the  matrix  equation 


jrfi!  =  x,  ‘/li 

|rf2i  o/Jk1 


where  [xj  and  1x^1  are  the  4  bv  4  electric  and  magnetic  susceptibility  tensors 
(matrices)  respectively.  Substitution  of  equation  (13)  back  into  equation  (7)  gives 


Mx  M2  /,  0||Xf  Oj/J  =  4n  s. 


[/,]  L0M>JL0hL^  c  Is: 


The  4  by  4  dielectric  and  permeability  matrices,  je|  and  |p| ,  are  related  to  the  electric 
and  magnetic  susceptibility  matrices  by  the  equations 

(el  =  U\  +  4 (  x J  (15) 


(pi  =  [/|  +  4k  |  x^l  (16) 

where  |/j  is  the  4  by  4  identity  matrix.  With  the  use  of  equations  (15)  and  (16),  matrix 
equation  (14)  can  be  written  in  the  form 


M,£  M-, 

[/.' 

_ 

M2  M,g 

h 

c 

1 

n  | 
Ha 

For  a  medium  having  electrical  conductive  properties,  the  electric  source  vector  |s, 
appearing  in  equation  (17)  obeys  Ohm's  law 


Is,  i— icy  0  f. ! 

{s2|  lO  O  f. 
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where  [a]  is  the  4  bv  4  electrical  conductivity  matrix.  \:ote,  in  the  absence  of  magnetic 
charge,  the  magnetic  source  vector  |s;)  may  be  set  equal  to  zero.  With  the  use  of 
equation  (18),  matrix  equation  (17)  can  be  written  in  the  form 

iM.e  +  f v  A1-,  f.  o 

1  '  -  ,  1  =  (19) 

I  M,  /VI,  q  f-,]  o 

where  [el  is  a  4  by  1  null  vector  and  where  bv  definition 

[Yl=  —  |o|.  (20) 

c 

Equation  (19)  is  the  matrix  representation  of  the  Maxwell  field  equations  for  a  linear, 
conductive  medium  of  infinite  extent.  If  the  medium  is  isotropic  and  homogeneous  as 
well,  the  4  bv  4  matrices  |e| ,  ||i|  ,  and  [o|  appearing  in  equations  (19)  and  (20)  have 
the  simplified  diagonal  form 


r 

e 

0 

0 

0] 

4 
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0 

0 

o 
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0 

01 
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o  o 

Im-I  = 
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,  |g|  = 
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oj 
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(21) 

0 
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o 

oi 

[0 

0 

0 

ej 

L° 

0 

0 

|T 

NO 

0 

0 

For  this  case,  all  of  the  nonzero  main-diagonal  elements  in  each  of  these  4  bv  4  matrices 
are  equal  to  a  scalar  constant,  which  is  independent  of  position  and  time.  Hence,  the 
dielectric  medium  being  considered  here  may  be  simply  described  by  the  three  scalar 
constants  e,  g,  o . 

2.3  MULTILAYER  DIELECTRIC  MEDIUM 

We  now  consider  the  geometry  shown  in  figure  1.  This  figure  depicts  a  multilayer 
optical  medium  consisting  of  N  layers  sandwiched  between  a  left  and  right  semi-infinite 
subregion  of  space.  Each  layer  represents  a  dielectric  subregion  of  space.  Adjacent  layers 
are  separated  by  parallel  planar  interface  boundaries  orthogonal  to  the  z-axis.  The 
multilayer  medium  is  illuminated  by  an  electromagnetic  wave  incident  from  the  left. 
This  gives  rise  to  both  a  reflected  and  a  transmitted  electromagnetic  wave.  The  principal 
objective  in  solving  a  problem  of  this  nature  is  to  determine  the  characteristics  of  the 
reflected  and  transmitted  waves  in  terms  of  the  properties  of  the  incident  wave  as  well  as 
the  optical  and  geometrical  properties  of  the  multilayer  medium. 
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Layer 

0  1  2  3  •  •  •  N  N+l 


Figure  1.  Multilayer  geometry. 


For  a  multilayer  optical  medium,  the  three  matrices  appearing  in  equation  (21)  are 
now  a  function  of  the  spatial  coordinate  2.  We  can  write  these  matrices  in  the  form 


,V  +  i 

[el  = 

n  =0 

N+l 

[Hi  =  ^[HW:> 

n  =  0 

N  +  1 

1®I  =  ^[o!„^(z). 

n  =0 


(22) 


(23) 


(24) 
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The  matrices  [e|n,  [ (j. |  and  [ajn  represent  the  dielectric,  permeability,  and 
conductivity  matrices  associated  with  subregion  n.  The  functions  4Jn(zj,  hereafter 
denoted  partition  functions,  are  used  in  partitioning  space  nto  the  various  subregions 
depicted  in  figure  1.  These  partition  functions  are  defined  in  terms  of  the  Heaviside  step 
function  (Bracewell,  1965).  The  reader  is  referred  to  appendix  A,  which  contains  a  more 
detailed  description  of  the  partition  functions  as  well  as  their  properties. 

The  electromagnetic  field  vector  appearing  in  equation  (19)  may  also  be  partitioned 
in  a  similar  manner.  That  is. 


where  [/,  f2]  T  is  the  electromagnetic  field  vector  associated  with  subregion  n.  The 
superscript  T  refers  to  transpose.  Next,  we  substitute  equations  (22)  through  (25)  back 
into  equation  (19).  Using  the  multiplicative  property  (see  appendix  A)  of  the  partition 
functions  H/n(z)  gives  the  following  matrix  representation  of  the  Maxwell  field 
equations  for  a  dielectric  multilayer  optical  medium 


N  +  1 


n  =  0  L 


MjE  +  fy  M2 

I 

U 

Vfz)  = 

-  - 

0 

M2  M tP 

a 

f_2_ 

J 

n  J 

0 

where  by  definition 


MjC  +  fy  M2 

1  =  +  M2 

1 - 

S 

to 

S 

L _ 

m2 

1  n  — 

(26) 


(27) 


Notice  the  similarity  between  the  matrix  form  (26)  of  the  Maxwell  field  equations  for  a 
multilayer  medium  and  the  matrix  form  (19)  for  an  infinite  medium  without  boundaries. 
Next,  we  will  show  that  the  matrix  representation  (26)  of  the  Maxwell  field  equations  for 
a  multilayer  medium  has  the  electromagnetic  boundary  conditions  built  into  the 
representation.  There  is  no  need  to  consider  the  use  of  the  Divergence  Theorem  and 
Stoke's  Theorem  in  establishing  the  boundary  conditions  of  the  electromagnetic  fields  at 
the  various  planar  interfaces. 
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With  the  use  of  the  product  rule  from  calculus  (Wylie,  1953),  equation  (26)  takes  on 
the  following  form 


N  +  1 


n  =  0 


M,e  +  iy  M, 
M2  M,|i 


2J„ 


Nt  1 

+X^u) 

n  -  0 


Cj£  c2! 

7.|  _  i 1,1 

c,  c{(t j 

-  -i  n  • 

fj  '  j  ..j 

where  from  appendix  A 


(28) 


(29) 


The  matrices  [ C,  ]  and  [C2]  appearing  in  equation  (28)  are  given  by 


0 

0 

0 

o’ 

’  0 

-1 

0 

o' 

0 

0 

0 

0 

+  1 

0 

0 

0 
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0 

0 

0 

-1 

and 

l  C2]  = 

0 

0 

0 

0 

0 

0 

+  1 

0 

i 

1 

L° 

0 

0 

0 

(30) 


With  the  use  of  equations  (A.5),  (A.6),  and  (A .7)  from  appendix  A,  we  can  rewrite 
equation  (28)  in  the  form 


N  +  1 

2>vz) 


n  =  0 


M,£  +  iy  M,] 

7,1 

1 - 

<: 

NJ 

1  n  1 

h 

(31) 


x«<~ 


p  =  o 


C,e  C2 
c2  C,p 


p  +  1 1 — ' 


p  +  I 


C]£  C2 

C2  c,p 


0 

.i/d. 


It  is  of  interest  to  noie  that  the  first  sum  appearing  in  equation  (31)  is  nothing  more  than 
a  mathematical  statement  of  the  fact  that  the  Maxwell  field  equations  must  be  satisfied 
for  each  dielectric  subregion  n,  (n  =  0, 1,  2, N  +  1 ) .  The  second  sum  appearing  in 
equation  (31)  contains  the  boundary  conditions  that  must  be  satisfied  at  each  planar 
interface  boundary  z  =  zp  ,  (p  =  0, 1, 2, N)  separating  the  various  subregions  of 
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space.  For  equation  (31)  to  be  satisfied  at  all  points  in  space,  the  following  two  sets  of 
matrix  equations  must  be  satisfied: 


Maxwell  field  equations  in  each  subregion  n 


A4,e  +  jy 

< 

-  ~t 

0 

m2 

M,g_ 

n 

k 

n 

0 

(32) 


Boundary  conditions  at  each  planar  interface  z  =  zp 


Ci* 

C21 

/,! 

C,e  C2 ' 

\u 

rj 

u 
_ 1 

p 

hi 

p 

f 

n 

to 

o 

(p  =  0,  1.2 . N). 


(33) 


Next,  we  will  consider  the  case  of  monochromatic  plane-wave  solutions.  Appendix  B 
contains  a  detailed  mathematical  discussion  on  monochromatic  plane-wave  solutions  of 
matrix  equation  (19),  or  equivalently  matrix  equation  (32),  and  their  resulting  properties. 
In  addition,  appendix  C  contains  a  detailed  mathematical  discussion  on  the  boundary 
condition  equation  (33)  for  the  case  of  monochromatic  plane-wave  solutions.  It  is  highly 
advisable  at  this  point  for  the  reader  to  become  familiar  with  the  concepts  and  notation 
presented  in  these  appendices  before  continuing  on. 

2.4  ELECTROMAGNETIC  PLANE-WAVE  FIELDS 

Suppose  the  dielectric  multilayer  optical  medium  depicted  in  figure  1  is  illuminated 
by  a  monochromatic  electromagnetic  plane-wave,  hereafter  referred  to  as  simply  the 
incident  wave.  The  incident  wave  will  give  rise  to  both  a  transmitted  and  a  reflected 
monochromatic  electromagnetic  plane-wave  as  shown  in  figure  1.  In  addition,  as  a  result 
of  the  incident  wave  excitation,  each  layer  within  the  multilayer  medium  will  support 
two  monochromatic  electromagnetic  plane-waves  propagating  in  opposite  directions 
(Fowles,  1968).  Therefore,  the  electromagnetic  field  vectors  appearing  in  equations  (32) 
and  (33)  must  be  of  the  form 


i/,1 

~/i 

+ 

h 

h 

R 

h 

(34) 
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where  each  of  the  field  vectors  on  the  right-hand  side  of  equation  (34)  has  a 
monochromatic  plane-wave  representation  given  by  equation  (B.2)  in  appendix  B, 
namely 


/l 

_ 

foi 

A 

R 

f()2 

'A 

/oi 

A 

L 

fo 12 

exp  {i[r\T[k\ R} 

R 


exp  {i[r]T  [k]  L} . 

L 


(35) 


(36) 


The  subscripts  R  and  L  appearing  in  equations  (34)  through  (36)  denote  wave 
propagation  to  the  right  and  left  respectively.  Substitution  of  equation  (34)  back  into 
equation  (33),  with  the  use  of  equations  (35)  and  (36),  gives  the  following  equation 


G]  O 
O  C, 


/oi 

/o2 


exp(  +  izpkz)  + 


701 

Ml 


exp(~izpkz ) 


(37) 


exp  (~izpkz) 

The  mathematical  details  le  ding  up  to  equation  (37)  are  presented  in  appendix  C.  The 
above  4  by  4  matrices  [G,]  and  [G,]  are  given  by  equation  (C-4). 

The  electric  field  vector  [/01]  and  magnetic  field  vector  [/02]  appearing  in  equation 
(37)  are  not  independent,  but  are  in  fact  related  to  one  another  through  equation  (B.20)  in 
appendix  B.  With  the  use  of  equation  (B.20),  we  can  eliminate  the  magnetic  field  vectors, 
for  example,  from  equation  (37).  Doing  so  allows  us  to  cast  equation  (37)  into  the 
following  matrix  form 


G,  O 

O  C, 


/oi 

fo2 


exp(+izpkz)  + 


S„<>)  SI3<:P> 

r  i/oi  i  R 

i 

I/oi  ]  K 

521  (zp)  S22(zp) 

[I/oiJ  L 

P 

P 

S;i  (zp>  S22  Up) 

P*  1 

E/°-  ■  ( 

(38) 
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where 


S,1  <-> 

S|2(z) 

1 

1 

_S2I(z) 

S22  (z)_ 

G.F.U) 

G2F}(z)R 


G,F2(z)  ] 
G2F2 (z) L  j 


and 


[F,  (z)  ]  s  [/]  exp  i  +  izkz)  and  [F2  (2)  ]  s  [/]  exp(-izkz) 


(39) 


(40) 


[R]  a-iji  [a,  +  a2|R  and  [L]  =~‘J^  la,  +  a,| / .  (41) 

The  matrix  [/]  in  equation  (40)  is  again  the  4  by  4  identity  matrix.  The  4  by  4  direction 
cosine  matrices  [aj  and  [a,]  are  defined  by  equation  (B.9)  in  appendix  B.  The 
quantity  e  is  the  complex  dielectric  constant,  also  defined  in  appendix  B.  Matrix 
equation  (38)  summarizes  the  relationship  between  the  electric  field  vectors  associated 
with  the  electromagnetic  waves  propagating  to  the  right  and  the  left  in  layers  p  and  p+1 
evaluated  at  the  planar  interface  z  =  zp. 

2.5  MULTILAYER  SYSTEM  MATRIX  OPERATOR 

We  will  next  generate  a  multilayer  system  matrix  operator  that  will  ultimately  allow 
us  to  determine  the  mathematical  relationship  between  the  incident,  reflected,  and 
transmitted  waves  depicted  in  figure  1.  First,  we  define  a  new  8  by  8  matrix  involving  the 
8  by  8  matrices  in  equation  (38),  namely 


- 

■ 

■ 

- 

-1 

r 

- 

w 

w 

__ 

s„ 

sn 

VV> 

S.2(ZP) 

.w 

K  (V 

p  +  l.p 

A1  ("V 

s,,  U„) 

P  +  1 

_^21  ^ ZP ^ 

S22(ZP) 

(42) 


The  superscript,  —  1,  denotes  the  multiplicative  matrix  inverse  operation.  Equation  (38) 
can  now  be  rewritten  in  the  form 


A  (z  J  A  (z  ) 
n  v  p  12  r 

R/0,1; 

A21(V  A22(2P] 

j  t/oi  1 , 

p  +  l.pw 
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Equation  (43)  allows  us  to  calculate  the  electric  field  vectors  in  subregion  p+1  in  terms  ot 
the  electric  field  vectors  in  subregion  p.  The  multilayer  system  matrix  operator  can  now 
be  computed  through  successive  matrix  multiplications  leading  to  the  following 
equation 


r 


12 


21 


*V> 

J  sv* 


n 

p  =  ,v 


A„  Up)  W 

A2i  (Zp)  A22(>} 


(44) 


p  t  !,p 


With  the  use  of  the  above  multilayer  system  matrix  operator,  we  can  now  express  the 
electric  field  vectors  in  the  right  semi-infinite  subregion  N+ 1  in  terms  of  the  electric  field 
vectors  in  left  semi-infinite  subregion  0.  The  result  is 


’[/oJki 

_  ! 

r 

ii 

r,2i 

r  t/o,  1 R 

N+  1 

/21 

b2j 

sys 

.[/oil, 

(45) 


For  simplicity  in  notation  we  let 

m  =  1  r/n,]R]0 

(46) 

l/,l  =  U/mltl„ 

(47) 

1/.)  = 

(48) 

The  three  vectors  in  equations  (46),  (47),  and  (48)  represent  the  incident,  reflected,  and 
transmitted  electric  field  vectors  respectively.  Based  on  the  illumination  scenario 
depicted  in  figure  1,  there  is  no  electromagnetic  wave  in  subregion  N+l  propagating  to 
the  left,  hence 

o 

II 

+ 

(49) 

With  the  use  of  equations  (46)  through  (49),  we  can  now  rewrite  equation  (45)  in  the  form 
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r  r 

11  12 

/, 

0 

r  r 

21  22  _ 

s  vs 

fr_ 

(50) 
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With  the  use  of  equation  (50),  we  can  now  solve  for  the  reflected  and  the  transmitted 
electric  field  vectors  in  terms  of  the  incident  electric  field  vector.  We  obtain  the  result 


[ U\  =  -l  r;‘  r21  ]  svs  i/,] 


[/,!  =  +  1  r  -  rp  r;’  i  f f,]. 


(51) 


(52) 


The  reflectance,  and  transmittance,  T,  coefficients  (Born  &  Wolf.  1965)  of  the 
multilayer  medium  can  now  be  computed  using  the  following  equations 


[/,]  H  [/,3 


(53) 


lf,]H[ft]  nt  cos  (6() 


T  : 


f/,-1  H  f  /,  1  ft  |  cos  (00 


(54) 


The  superscript,  h,  denotes  the  Hermitian  conjugate  operation.  The  quantities  ni  and  nt 
correspond  to  the  ordinary  indices  of  refraction  for  subregions  0  and  N  +  1  respectively. 
The  angles  0(  and  0(  represent  the  angle  of  incidence  and  angle  of  refraction  for  the 
subregions  0  and  N  +  1  respectively.  These  angles  and  indices  of  refraction  obey  the 
Snell's  law  relation 

tt.sin  (0.)  =  n,sin  (0().  (55) 

Equation  (54)  tacitly  assumes  that  the  subregions  0  and  N+l  are  nonconducting. 
Typically  these  two  subregions  correspond  to  either  vacuum,  air,  or  glass.  There  are  no 
restrictions  placed  on  the  optical  properties  of  the  dielectric  layers  constituting  the 
multilayer  medium.  Equations  (53)  and  (54)  are  the  quintessential  equations  necessary 
for  calculating  the  reflection  and  transmission  characteristics  of  the  multilayer  medium. 
It  is  noted  that  the  wavelength,  angle  of  incidence,  and  polarization  state  of  the  incident 
wave  can  be  arbitrarily  chosen.  These  equations  serve  as  the  basis  for  the  applications 
considered  in  the  following  section. 
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3.0  APPLICATIONS 


Based  on  the  8  by  8  matrix  representation  of  electromagnetic  theory  being  developed 
by  the  author,  a  computer  software  package  is  currently  being  written  and  tested  at 
NRaD  to  handle  a  variety  of  scenarios  involving  electromagnetic  wave  propagation 
through  matter.  This  software  package  is  referred  to  as  the  MATURE  Program.  MATURE 
is  the  acronym  for  Matrix  Approach  To  Understanding  Relativistic  Electrodynamics. 
MATURE  is  written  in  MATLAB  code  for  a  Sun  4  SPARCstation  2  workstation.  One 
feature  of  the  MATURE  Program  is  its  ability  to  solve  electromagnetic  wave  propagation 
problems  involving  multilayer  optical  structures.  The  required  inputs  to  the  MATURE 
Program  for  problems  of  this  nature  are  based  on  the  following  statements: 

1.  Enter  total  number  of  optical  layers,  N. 

2.  Enter  thickness,  T  ,  of  each  layer  in  nanometers  (nm)  for  p  =  7,  2,...,  N. 

3.  Enter  ordinary  refractive  index,  np,  of  each  subregion  p  =  0,1,  2,...,  N,  N+l. 

4.  Enter  attenuation  index,  k  of  each  subregion  p  =  0,1,  2,...,  N,  N+l. 

5.  Enter  angle  of  incidence,  9t,  of  incident  plane-wave  in  degrees  (deg). 

6.  Enter  wavelength,  X,  of  incident  plane-wave  in  nm. 

7.  Specify  polarization  state,  Transverse  Electric  (TE)  or  Transverse  Magnetic  (TM), 
of  incident  plane- wave. 

8.  Determine  the  variable  of  interest,  namely  wavelength,  X,  angle  of  incidence,  0., 
or  thickness,  Tn,  of  the  nth  layer. 

9.  If  the  variable  of  interest  is  the  wavelength,  then 

a.  Specify  minimum  value  of  wavelength,  Xmin,  in  nm, 

b.  Specify  maximum  value  of  wavelength,  Xmax,  in  nm, 

c.  Specify  incremental  value  of  wavelength,  A^.,  in  nm. 

10.  If  the  variable  of  interest  is  the  angle  of  incidence,  then 

a.  Specify  minimum  value  of  angle  of  incidence,  0|  mm,  in  deg, 

b.  Specify  maximum  value  of  angle  of  incidence,  0(  max,  in  deg, 

c.  Specify  incremental  value  of  angle  of  incidence,  A0;,  in  deg. 

11.  If  the  variable  of  interest  is  the  thickness  of  the  nth  layer,  then 

a.  Specify  layer  number,  n,  of  interest, 

b.  Specify  minimum  thickness,  Tn  min,  of  nth  layer  in  nm, 

c.  Specify  maximum  thickness,  Tn  max,  of  nth  layer  in  nm, 

d.  Specify  incremental  thickness,  ATn,  of  nth  layer  in  nm. 

The  outputs  to  the  MATURE  Program  are  the  reflectance,  and  transmittance,  T,  as  a 
function  of  the  variable  of  interest. 
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With  the  use  of  the  MATURE  Program,  a  number  of  examples  are  considered  in  this 
section  that  illustrate  the  utility  of  the  8  by  8  matrix  representation  in  solving 
electromagnetic  wave  propagation  problems  involving  multilayer  optical  media.  The 
examples  considered  involve:  (1)  reflection  and  refraction  at  a  single  interface  separating 
two  nonconducting  dielectric  media,  (2)  evanescent  wave  coupling  from  one  dielectric 
medium  into  another  via  a  vacuous  thin  film,  (3)  the  reflection  and  transmission 
properties  of  a  single  conducting  metal  thin  film,  (4)  the  reflection  and  transmission 
characteristics  of  a  multilayer  antireflection  coating  structure,  (5)  the  reflection  and 
transmission  properties  of  a  multilayer  high-reflectance  mirror,  and  (6)  the  reflection  and 
transmission  characteristics  of  a  multilayer  interference  band-pass  filter. 

3.1  SINGLE  SURFACE  REFLECTION  AND  REFRACTION 

For  our  first  example  we  consider  reflection  and  refraction  at  a  planar  interface 
separating  two  nonconducting  dielectric  media.  This  corresponds  to  a  geometry  in 
which  N  =  0  in  figure  1.  It  is  convenient  to  consider  two  special  cases.  The  first  case  is  that 
in  which  the  electric  field  vector  of  the  incident  wave  is  parallel  to  the  planar  interface. 
This  is  commonly  referred  to  as  the  transverse  electric  or  TE  polarization  case.  The 
second  case  is  that  in  which  the  magnetic  field  vector  of  the  incident  wave  is  parallel  to 
the  interface.  This  is  known  as  the  transverse  magnetic  or  TM  polarization  case.  The 
general  case  is  handled  by  using  appropriate  linear  combinations  of  these  two  special 
cases.  Figure  2  shows  the  results  obtained  from  the  MATURE  Program  for  a  rare-to- 
dense  (external  reflection)  scenario  (Fowles,  1968).  The  left  semi-infinite  subregion  of 
space  in  figure  1  is  a  vacuum  and  the  right  semi-infinite  subregion  is  a  nonconducting 
medium  with  ordinary  refractive  index  equal  to  1.5.  The  wavelength  of  the  incident 
wave  was  500  nm  and  the  angle  of  incidence  was  chosen  as  the  variable  of  interest.  It  is 
apparent  from  figure  2  that,  for  the  TM  case,  the  reflectance  is  zero  at  the  Brewster  angle 
(polarization  angle)  value  of  56.3°  predicted  by 

r«,l  r  1.5  n 

0fc  =  atan  —  =  atan  —  .  (56) 

L  wo  j  l  IQ. 

0°  is  equal  to  4  percent,  which  is  in  agreement  with 

(57) 

used  to  predict  (Fowles,  1968)  reflectance  at  normal  incidence.  Figure  3  illustrates  the 
results  obtained  for  a  dense-to-rare  (internal  reflection)  scenario  (Fowles,  1968).  For  this 


Also  note  that  the  reflectance  at  0  = 

I 

the  formula 
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Figure  2.  Single  surface  (rare-to-dense). 
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Figure  3.  Single  surface  (dense-to-rare). 


scenario,  the  right  semi-infinite  subregion  of  space  is  a  vacuum  and  the  left  semi-infinite 
subregion  is  a  nonconducting  medium  with  ordinary  retractive  index  equal  to  15 
Again,  the  wavelength  was  set  equal  500  nm  in  the  MATURE  Program  and  the  angle  of 
incidence  was  chosen  as  the  variable  of  interest.  The  Brewster  angle  for  this  scenario  is 
equal  to  33.7°.  The  reflectance  at  normal  incidence  is  again  4  percent.  In  addition,  we 
observe  that  total  internal  reflectance  occurs  for  angles  of  incidence  greater  than  the 
critical  angle  (Fowles,  1968)  of  41.8°  predicted  by  the  formula 


'«T 

.  n.oi 

=  asin 

_”o_ 

=  asm  | 

!  l1-3_ 

(58) 


The  reflectance  and  transmittance  plots  in  both  figures  2  and  3  are  in  excellent  agreement 
with  results  found  in  standard  textbooks  on  optics,  e  g..  Born  &  Wolf  (1965)  and  Fowles 
(1968),  treating  this  subject. 


3.2  EVANESCENT  WAVE  COUPLING 


In  the  next  example,  we  consider  the  transport  of  electromagnetic  energy  between 
two  semi-infinite  nonconducting  dielectric  media  separated  by  a  vacuous  thin  film.  This 
corresponds  to  N  =  1  in  the  geometry  of  figure  1.  For  this  example,  the  ordinary 
refractive  index  of  each  dielectric  medium  was  set  equal  to  1.5.  The  wavelength  of  the 
incident  wave  was  set  equal  to  500  nm.  Only  the  TF  polarization  case  was  considered. 
The  thickness  of  the  vacuous  thin  film  was  varied  between  0  and  1000  nm.  Three 
different  test  cases  were  considered  corresponding  to  angles  of  incidence  of  0°,  38°,  and 
45°,  respectively.  Figure  4  summarizes  the  results  obtained  using  the  MATURE  Program. 
Note  that  the  critical  angle  for  this  example  is  equal  to  41.8°.  For  angles  of  incidence  less 
than  the  critical  angle,  the  vacuous  thin  film  behaves  as  an  optical  cavity.  The  periodic 
variations  in  both  the  reflectance  and  transmittance  of  the  multilayer  optical  structure, 
with  thin-film  thickness,  is  clearly  obvious  from  the  plots  in  figure  4.  However,  for  angles 
of  incidence  exceeding  the  critical  angle,  the  periodic  behavior  of  the  reflectance  and 
transmittance  of  the  multilayer  structure  ceases.  This  periodic  behavior  has  been 
replaced  by  an  exponential  decrease  (increase)  in  the  transmittance  (reflectance)  of  the 
multilayer  structure  with  increasing  values  of  the  thin-film  thickness.  Notice  how  the 
optical  structure  becomes  totally  reflecting  for  all  practical  purposes  beyond  a  thin-film 
thickness  of  about  600  nm  when  the  angle  of  incidence  is  equal  to  45°.  The  behavior  of  a 
structure  of  this  nature  for  angles  of  incidence  exceeding  the  critical  angle  is  attributed  to 
a  phenomena  called  evanescent  wave  coupling.  For  a  detailed  discussion  on  evanescent 
wave  coupling,  the  reader  is  referred  to  Fowles  (1968). 
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Figure  4.  Vacuous  thin  film  (TE  polarization). 
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3.3  SINGLE  CONDUCTING  THIN  FILM 


For  our  next  example,  we  consider  a  multilayer  optical  medium  consisting  of  a  single 
conducting  metallic  thin  film  separating  two  semi-infinite  dielectric  subregions  of  space. 
This  corresponds  to  N  =  1  in  the  geometry  of  figure  1 .  There  are  two  scenarios  we  wish  to 
consider.  In  the  first  scenario,  the  left  semi-infinite  subregion  of  space  in  figure  1  is  a 
vacuum  and  the  right  semi-infinite  subregion  is  a  dielectric  medium  with  ordinary 
refractive  index  equal  to  1.5.  The  ordinary  refractive  index  of  the  metal  film  is  3.5.  The 
attenuation  index  of  the  metallic  thin  film  was  set  equal  to  three  different  values,  namely 
0.0,  0.1,  and  0.2.  The  thin-film  structure  was  normally  illuminated  by  a  plane-wave  of 
wavelength  500  nm.  For  normal  incidence  there  is  no  distinction  between  the  TE  and  TM 
polarization  cases.  The  variable  of  interest  was  the  thin-film  thickness.  The  thickness  was 
varied  between  0  and  300  nm.  In  figures  5,  6,  and  7  are  plots  of  the  reflectance  and 
transmittance  of  the  metallic  thin-film  structure  as  a  function  of  thin-film  thickness  for 
the  attenuation  index  values  of  0.0,  0.1,  and  0.2  respectively.  The  results  summarized  in 
these  three  figures  are  in  excellent  agreement  with  similar  results  reported  by  Born  and 
Wolf  (1965). 

The  second  scenario  considered  involves  a  metallic  thin  film  separating  two  vacuous 
semi-infinite  subregions  of  space.  The  ordinary  refractive  and  attenuation  index  values, 
(«,,  Kj ) ,  chosen  for  the  thin  film  correspond  to  values  for  real  metals.  In  particular,  the 
following  metals  were  considered:  silver  (0.200,  17.200);  aluminum  (1.440,  3.632);  gold 
(0.470,  6.021);  copper  (0.620,  4.145);  and  iron  (1.510,  1.079).  The  thin-film  structure  was 
illuminated  at  normal  incidence  by  a  plane-wave  of  wavelength  589.3  nm.  The  variable 
of  interest  was  the  thickness  of  the  thin  film  that  was  varied  between  0  and  100  nm. 
Results  obtained  for  the  five  metallic  thin  films  considered  are  shown  in  figures  8 
through  12.  Each  of  these  figures  clearly  indicates  that  the  transmittance  of  a  metallic 
thin  film  approaches  zero  as  the  thin-film  thickness  increases.  In  addition,  the  reflectance 
of  a  metallic  thin  film  approaches  a  maximum  value  with  increasing  thin-film  thickness. 
The  maximum  value  of  the  reflectance  is  predicted  (Born  &  Wolf,  1965)  by  the  equation 


K 


ax 


n]  ( 1  +  ic~)  +  1  -  2 h, 
tt]  ( 1  +  xq)  +  1  +  2 


(59) 


The  maximum  values  of  the  reflectance  obtained  for  the  various  thin  films  considered 
are  from  equation  (59):  silver  (0.94);  aluminum  (0.83);  gold  (0.82):  copper  (0.73);  and  iron 
(0.33). 
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Figure  5.  Metallic  thin  film,  attenuation  index  value  of  0.0. 
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Figure  6.  Metallic  thin  film,  attenuation  index  value  of  0.1. 
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Figure  7.  Metallic  thin  film,  attenuation  index  value  of  0.2. 
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Figure  8.  Silver  tifn  film. 
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Figure  9.  Aluminum  thin  film. 
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Figure  10.  Gold  thin  film. 
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Figure  1 1 .  Copper  thin  film. 


28 


REGION 

THICKNESS 

(nm) 

REFRACTIVE 

INDEX 

ATTENUATION 

INDEX 

LEFT 

1.000 

0.000 

1 

100.0 

1.510 

1.079 

RIGHT 

1.000 

0.000 

Incidence  Angle:  0.0  (deg)  Wavelength:  589.3  (nm) 


REFLECTANCE  (%) 


Thickness  (nm) 


Figure  12.  Iron  thin  film. 


3.4  ANTIREFLECTION  COATINGS 

When  an  electromagnetic  wave  impinges  upon  an  optical  surface  of  a  prism  or  lens,  a 
portion  of  the  wave  is  reflected  at  the  surface.  This  reflection  loss  can  have  deleterious 
effects  on  the  overall  transmission  characteristics  of  the  optical  device.  Reflection  losses 
reduce  the  overall  transmission  of  optical  devices  and  instruments  and  often  produce 
unwanted  ghost  images.  The  use  of  antireflection  coatings  on  optical  surfaces  is  a 
practical  method  for  reducing  these  undesirable  effects.  A  detailed  treatment  on  the 
theory  of  antireflection  coatings  may  be  found,  for  example,  in  Baumeister  u  5), 
Musset  and  Thelen  (1970),  and  MacLeod  (1986). 

We  now  consider  a  multilayer  optical  structure  consisting  of  anywhere  from  one  to 
four  thin-film  layers  for  the  purpose  of  reducing  the  reflection  loss  at  the  interface 
between  two  dielectric  media.  In  particular,  we  will  consider  the  case  when  the  left  semi¬ 
infinite  subregion  of  space  in  figure  1  is  vacuum  and  the  right  semi-infinite  subregion  is 
glass  with  an  ordinary  refractive  index  equal  to  1.52.  The  multilayer  optical  structure  is 
illuminated  by  a  plane-wave  at  normal  incidence.  The  wavelength  was  chosen  as  the 
variable  of  interest  with  values  ranging  over  visible  portion  of  the  electromagnetic 
spectrum  from  400  nm  (blue)  to  700  nm  (red).  Figures  13  through  16  show  the  reflectance 
and  transmittance  plots  obtained  when  antireflection  coatings  are  employed.  The 
thicknesses  and  optical  properties  of  the  various  layers  comprising  the  various 
multilayer  structures  considered  are  summarized  in  these  figures.  It  is  obvious  from 
these  figures  that  the  reflection  loss  at  a  single  surface  for  normal  incidence  is 
approximately  4  percent  when  the  surface  is  uncoated.  Using  a  single  antireflection 
coating  (figure  13)  can  reduce  the  losses  to  approximately  2  percent  over  the  visible 
spectrum.  Two  coatings  (figure  14)  can  reduce  the  reflection  losses  to  approximately  1 
percent.  Three  and  four  coatings  (figures  15  and  16)  can  further  reduce  and  even 
eliminate  these  reflection  losses  over  a  substantial  portion  of  the  visible  spectrum.  The 
MATURE  Program  computer  results  presented  in  figures  13  through  16  are  in  excellent 
agreement  with  results  obtained  by  Musset  and  Thelen  (1970)  for  the  same  multilayer 
structures. 

Antireflection  coatings  can,  of  course,  be  used  to  eliminate  unwanted  reflections  in 
other  portions  of  the  electromagnetic  spectrum  as  well.  Figure  17  summarizes  the  results 
obtained  from  the  MATLAB  4.0  computer  program  for  wavelengths  in  the  infrared 
portion  of  the  spectrum  for  a  two-layer  antireflection  coating.  These  results  are  also  in 
excellent  agreement  with  those  reported  by  MacLeod  (1986). 


30 


REGION 

THICKNESS 

(nm) 

REFRACTIVE 

INDEX 

ATTENUATION 

INDEX 

LEFT 

1.000 

0.000 

1 

92.4 

1.380 

0.000 

RIGHT 

1.520 

0.000 

Incidence  Angle:  0.0  (deg) 


REFLECTANCE  (%)  TRANSMITTANCE  (%) 


Wavelength  (nm)  Wavelength  (nm) 

Figure  13.  Antireflection  coating  (1  layer). 
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Figure  14.  Antireflection  coating  (2  layers). 
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Figure  15.  Antireflection  coating  (3  layers). 
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Figure  16.  Antireflection  coating  (4  layers). 
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Figure  17.  Antireflection  coating  (IR). 


3.5  HIGH-REFLECTANCE  COATINGS 


Overcoating  a  glass  substrate  with  a  multilayer  dielectric  thin-film  coating  is  one 
practical  method  of  fabricating  high-reflectance  mirrors  used  in  both  Fabry-Perot 
interferometers  and  laser  cavities.  Often,  these  optical  devices  require  mirrors  that  not 
only  have  a  reflectance  close  to  100  percent,  but  which  also  transmit  a  substantial  portion 
of  the  electromagnetic  energy  that  is  not  reflected.  Metallic  thin  films  of  silver, 
aluminum,  or  gold  seldom  attain  a  reflectance  that  exceeds  90  percent  in  the  visible 
portion  of  the  electromagnetic  spectrum.  However,  multilayer  dielectric  mirrors  used  for 
this  purpose  can  attain  a  reflectance  of  over  99  percent  (Baumeister,  1965). 

One  type  of  multilayer  dielectric  coating  used  in  the  fabrication  of  high-reflectance 
mirrors  is  the  quarter-wave  stack  (Baumeister,  1965).  Basically,  the  quarter-wave  stack  is 
a  multilayer  structure  where  each  layer  has  the  same  optical  thickness,  and  with 
refractive  indices  that  alternate  between  two  values.  The  optical  thickness,  not  to  be 
confused  with  geometrical  thickness  of  a  layer,  is  simply  the  product  of  the  refractive 
index  of  the  layer  with  its  geometrical  thickness. 

We  now  consider  a  multilayer  structure  having  a  quarter-wave  stack  geometry  for 
the  purpose  of  increasing  the  reflectance  at  the  interface  between  two  dielectric  media. 
Four  cases  are  considered:  (1)  single-layer;  (2)  three-layer;  (3)  five-layer;  and  (4)  a  15- 
layer  stack.  We  will  consider  the  case  when  the  left  semi-infinite  subregion  of  space  in 
figure  1  is  vacuum  and  the  right  semi-infinite  subregion  is  glass  with  an  ordinary 
refractive  index  equal  to  1.51.  Let  the  multilayer  quarter-wave  stack  be  illuminated  by  a 
plane-wave  at  normal  incidence.  The  wavelength  was  chosen  as  the  variable  of  interest. 
The  range  of  wavelength  values  varied  between  450  nm  and  650  nm.  Figures  18  through 
21  show  the  reflectance  and  transmittance  plots  obtained  when  high-reflectance  coatings 
are  employed.  The  geometrical  thicknesses  and  optical  properties  of  the  various  layers 
comprising  the  quarter-wave  stack  considered  are  summarized  in  these  figures.  It  is 
obvious  from  these  figures  that  the  reflectance  at  a  single  surface  for  normal  incidence  is 
approximately  4  percent  when  the  surface  is  uncoated.  Using  a  single  high-reflectance 
coating  (figure  18)  can  increase  the  reflectance  to  over  30  percent  over  the  spectral  range 
of  interest.  Three  coatings  (figure  19)  can  increase  the  reflectance  to  over  50  percent  and 
five  coatings  (figure  20)  can  increase  the  reflectance  to  over  70  percent.  A  15  layer 
quarter-wave  stack  (figure  21)  can  increase  the  reflectance  to  over  99  percent  over  a 
substantial  portion  of  the  visible  spectrum.  The  results  of  the  MATURE  Program 
presented  in  figures  18  through  21  are  in  excellent  agreement  with  results  reported  by 
Musset  and  Thelen  (1970)  for  the  same  quarter-wave  stack  structures. 
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Figure  18.  High-reflectance  coating  (1  layer). 


37 


REGION 

THICKNESS 

(nm) 

REFRACTIVE 

INDEX 

ATTENUATION 

INDEX 

LEFT 

1.000 

0.000 

1 

57.8 

2.300 

0.000 

2 

96.4 

1.380 

0.000 

3 

57.8 

2.300 

0.000 

RIGHT 

1.510 

0.000 

Incidence  Angle:  0.0  (deg) 


REFLECTANCE  (%)  TRANSMITTANCE  (%) 


Wavelength  (nm)  Wavelength  (nm) 

Figure  19.  High-reflectance  coating  (3  layers). 
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Figure  20.  High-reflectance  coating  (5  layers). 
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Figure  21.  High-reflectance  coating  (15  layers). 
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3.6  MULTILAYER  INTERFERENCE  BAND-PASS  FILTER 

Multilayer  interference  band-pass  filters  are  often  used  in  place  of  a  prism  or 
diffraction  grating  in  spectrophotometric  systems.  These  filters,  for  a  given  resolving 
power,  are  the  most  luminous  instrument  available.  The  advantage  of  using  only 
dielectric  films,  as  opposed  to  metallic  films,  is  that  the  peak  transmittance  is  high  and  a 
narrowband  pass  can  be  achieved.  A  detailed  discussion  on  the  fabrication  and  use  of 
multilayer  interference  filters  can  be  found  in  Hernandez  (1986). 

The  quarter-wave  stack  previously  mentioned  in  our  discussion  on  high-reflectance 
coatings  also  finds  use  in  the  fabrication  of  multilayer  interference  band-pass  filters.  One 
geometry  involves  sandwiching  a  half-wave  spacer  between  two  identical  quarter-wave 
stacks  (Baumeister,  1965).  We  now  consider  a  multilayer  structure  having  this  geometry. 
Three  cases  are  considered:  (1)  a  7-layer;  (2)  an  11-layer;  and  (3)  a  15-layer  filter.  We 
consider  the  scenario  for  which  the  left  semi-infinite  subregion  of  space  in  figure  1  is 
vacuum  and  the  right  semi-infinite  subregion  is  glass  with  an  ordinary  refractive  index 
equal  to  1.51.  Let  the  multilayer  filter  be  illuminated  by  a  plane- wave  whose  wavelength 
is  varied  between  500  nm  and  600  nm.  Shown  in  figures  22  through  24  are  reflectance 
and  transmittance  plots  obtained  for  the  7-layer,  11 -layer,  and  15-layer  filters  respectively 
for  the  case  of  normal  incidence.  The  geometrical  thicknesses  and  optical  properties  of 
the  layers  comprising  these  various  multilayer  filters  are  also  summarized  in  these 
figures.  It  is  obvious  from  these  figures  that  the  band  pass  of  an  interference  filter 
decreases  as  the  number  of  layers  in  each  of  the  quarter-wave  stacks  increases.  Shown 
next  in  figures  25  and  26  is  the  effect  angle  of  incidence  has  on  the  reflectance  and 
transmittance  characteristics  of  an  11-layer  interference  filter.  Plots  for  angles  of 
incidence  equal  to  0°,  30°,  and  45°  are  shown.  The  results  presented  in  figure  25  are  for 
an  incident  wave  having  a  TE  polarization  state,  whereas  those  in  figure  26  correspond 
to  an  incident  wave  having  a  TM  polarization  state.  It  is  obvious  from  these  figures  how 
the  band  pass  of  the  interference  filter  is  shifted  to  shorter  wavelengths  as  the  angle  of 
incidence  is  increased.  We  see  that  the  interference  filter  characteristics  also  depend  on 
the  polarization  state  of  the  incident  plane-wave. 

Shown  in  figure  27  are  reflectance  and  transmittance  plots  obtained  using  the 
MATURE  Program  of  a  three-laver  Fabry-Perot  etalon  where  two  identical  silver  thin 
films  are  used  in  conjunction  with  a  single  dielectric  spacer.  Notice,  as  a  result  of 
absorption  within  the  metallic  films  of  the  etalon,  how  the  overall  transmittance  of  the 
etalon  is  greatly  reduced  in  comparison  to  that  of  an  all-dielectric  interference  filter.  This 
demonstrates  the  advantage  of  using  only  dielectric  thin  films  in  interference  filters. 
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Figure  22.  Band-pass  filter  (7  layers). 
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Figure  23.  Band-pass  filter  (1 1  layers). 
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Figure  24.  Band-pass  filter  (1 5  layers). 
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Figure  25.  Band-pass  filter  (1 1  layers,  TE  polarization). 
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Figure  26.  Band-pass  filter  (1 1  layers,  TM  polarization). 
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Figure  27.  Fabry-Perot  etalon  with  silver  mirrors. 
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4.0  SUMMARY  AND  CONCLUSIONS 


Presented  in  this  publication  is  a  new  approach,  based  entirely  on  the  use  of  matrices, 
for  solving  electromagnetic  wave  propagation  problems  involving  multilaver  dielectric 
media  with  planar  boundaries.  This  new  approach  is  based  upon  a  reformulation  of  the 
traditional  vector  form  of  the  Maxwell  field  equations  into  a  single  matrix  equation 
containing  square  8  by  8  differential  matrix  operators.  Boundary  conditions  are  taken 
into  account  by  partitioning  both  the  electromagnetic  wavefunctions  and  the  dielectric, 
permeability,  and  conductivity  matrices  describing  the  optical  medium  using 
combinations  of  Heaviside  step  functions  called  partition  functions.  Using 
monochromatic  plane- wave  representations  of  the  electromagnetic  field,  we  were  able  to 
develop  an  8  by  8  system  operator  description  of  wave  propagation  through  the  entire 
multilayer  medium.  This  operator  description  led  to  a  set  of  equations  for  predicting 
both  the  reflection  and  transmission  characteristics  of  a  multilayer  medium  with  planar 
boundaries  illuminated  by  a  monochromatic  plane-wave  electromagnetic  radiation. 

Because  of  the  power  of  matrix  operations,  many  of  the  fundamental  relationships 
and  results  in  electromagnetic  wave  propagation  theory  are  easily  derived  without  the 
need  for  the  usual  plethora  of  vector  calculus  identities  that  have  become  the  standard  in 
these  derivations:  Stokes'  theorem,  the  Divergence  theorem,  the  formula  for  the  curl  of  a 
curl,  etc.  Instead,  mathematical  operations  involving  matrix  multiplication,  matrix 
inversion,  and  eigensolutions  of  matrices  can  be  employed.  Readily  available  off-the- 
shelf  computer  software  packages,  like  MATLAB  for  numerical  computations  and 
MATHEMAT1CA  for  symbolic  manipulations,  are  well  suited  for  these  matrix 
operations.  A  computer  software  package  entitled  MATURE  is  currently  being  written 
and  tested  at  NRaD  to  handle  a  variety  of  scenarios  involving  electromagnetic  wave 
propagation.  MATURE  is  the  acronym  for  Matrix  Approach  To  Understanding 
Relativistic  Electrodynamics.  The  MATURE  computer  software  program  is  being  written 
using  MATLAB  code  for  use  on  a  Sun  4  SPARCstation  2  workstation. 

With  the  present  version  of  the  MATURE  Program,  several  numerical  examples  were 
considered  that  illustrate  the  power  and  utility  of  using  the  8  by  8  matrix  representation 
in  solving  wave  propagation  problems.  Examples  considered  included:  (1)  reflection  and 
refraction  at  a  single  interface  separating  two  dielectric  media,  (2)  evanescent  wave 
coupling  from  one  dielectric  medium  into  another  via  a  vacuous  thin  film,  (3)  the 
reflection  and  transmission  properties  of  a  single  conducting  thin  metallic  film,  (4)  the 
reflection  and  transmission  characteristics  of  a  multilayer  antireflection  coating 
structure,  (5)  the  reflection  and  transmission  properties  of  a  multilayer  high-reflectance 
mirror,  and  (6)  the  reflection  and  transmission  characteristics  of  a  multilayer  interference 
band-pass  filter.  The  MATURE  Program  is  quite  general  in  that  the  reflection  and 
transmission  characteristics  of  the  multilayer  medium  can  be  determined  for  arbitrarily 
chosen  wavelengths,  angles  of  incidence,  and  polarization  states  of  the  incident 
electromagnetic  radiation. 
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APPENDIX  A 


PARTITION  FUNCTIONS 

The  partition  functions,  first  introduced  in  section  2.3  of  this  report,  are  defined  in 
this  appendix.  Recall  that  the  multilayer  optical  medium  depicted  in  figure  1  consists  of 
a  finite  number,  N,  of  planar  layers  each  having  a  finite  thickness  as  measured  along  the 
z-direction.  All  layer  interfaces  are  parallel  as  well  as  being  orthogonal  to  the  2-axis.  Each 
of  these  layers  is  labelled  by  an  integer,  n,  where  n  =  1,2,  ...,N.  To  the  left  of  the 
multilayer  medium  is  a  semi-infinite  subregion  of  space  denoted  by  the  label  0.  And  to 
the  right  of  the  multilayer  medium  is  a  second  semi-infinite  subregion  of  space  denoted 
by  the  label  N+l.  Associated  with  each  of  the  layers  comprising  the  multilayer  medium 
and  the  two  semi-infinite  subregions  of  space  is  a  partition  function.  These  partition 
functions  are  defined  by  equations  (A.l)  through  (A.3)  below.  Those  properties  of  the 
partition  functions  pertinent  to  this  investigation  are  contained  in  equations  (A.4) 
through  (A. 8).  It  is  noted  that  the  partition  functions  are  simply  Heaviside  unit  step 
functions  (Bracewell,  1965)  or  linear  combinations  thereof. 

The  partition  function  associated  with  the  semi-infinite  subregion  of  space  labelled  0 
is  defined  by  the  following  equation 

r  1  2  <  20 

¥0(z)  =  |  (A.l) 

^  0  otherwise. 

Associated  with  each  layer,  n,  of  the  multilayer  optical  medium  is  a  partition  function 
defined  by 


^(2)=  (n  -  1,  2, 3, ...,  N)  (A. 2) 

^  0  otherwise 

The  semi-infinite  subregion  of  space  labelled  N+l  is  represented  by  a  partition  function 
defined  by 


0  otherwise. 


(A.3) 


A-l 


All  of  the  partition  functions  defined  above  satisfy  the  following  property 


(A.4) 

where  8mn  is  the  Kronecker  delta  (Bracewell,  1965).  Each  partition  function  4>n  (z)  has  a 
first  derivative  with  respect  to  the  variable  z,  denoted  by  the  function  <J>n  (z)  . 
Mathematically 

<Dn(z)  H-fy  (2)  (n  =  0, 1,  2, 3, ....  N,  N  +  1) .  (A. 5) 

"  dz  " 

Using  the  fact  that  the  first  derivative  of  the  Heaviside  unit  step  is  the  Dirac  delta 
function  (Bracewell,  1965),  it  can  be  easily  shown  that 

<D0(z)  =  -8(z-z0)  (A.6) 

<*>„ (z)  =  +5(z-z„_1) -5(z-z„)  (rt  =  1, 2, 3, ....  N)  (A.7) 

°N+l(2)  =  +5(2-zn)  (A.8) 

where  5(z)  represents  the  Dirac  delta  function. 
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APPENDIX  B 


ELECTROMAGNETIC  PLANE-WAVE  SOLUTIONS 


Monochromatic  electromagnetic  plane-wave  solutions  satisfying  the  Maxwell  field 
equations  are  of  vital  importance  to  our  understanding  of  wave  propagation  through  a 
multilayer  optical  medium.  We  start  with  the  matrix  representation  of  the  Maxwell  field 
equations  for  an  unbounded  dielectric  medium  as  summarized  bv  equation  (19)  and 
rewritten  here  in  equation  (B.l)  for  convenience. 


M1e  +  iy  M2 

>1 

0 

CN 

_ 1 

h. 

0 

For  monochromatic  electromagnetic  plane-wave  fields,  the  electromagnetic  field  vector 
appearing  in  equation  (B.l)  can  be  expressed  in  the  form 


u 

h 


/oi 

/o2 


exp{i[r]T[k)}. 


(B.2) 


The  4  by  1  vectors  [/01]  and  [f^]  are  both  constant  and  describe  the  polarization 
properties  of  the  plane-wave  field.  With  the  use  of  equation  (10),  they  can  be  expressed 
in  the  form 
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(B.3) 


The  relativistic  4-vectors  f  r]  and  [A:]  appearing  in  equation  (B.2)  are  defined  by 
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(B.4) 
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where 


x  =  ict  and  k=  — .  (B.5) 

c 

Again,  the  imaginary  unit  i  is  equal  to  J-i  and  c  represents  the  speed  of  light  in 
vacuum.  The  quantity  to  represents  the  angular  frequency  of  the  monochromatic  plane- 
wave.  The  superscript,  T,  appearing  in  equation  (B.2)  denotes  transpose.  The  dot  (scalar) 
product  of  the  vectors  [  r]  and  [fc]  is  given  by 

[r]T  [k]  =  xkx  +  ykv  +  zkz-  at  =  k»r-at.  (B.6) 

The  three-dimensional  wave  vector  k  above  can  be  expressed  in  the  form 

k  =  (kx,ky,kz)  =  k(ax,ay,az)  =  ka  (B.7) 


where  the  scalar  k  is  the  wavenumber.  The  unit  vector  a  has  as  its  components  the 
direction  cosines  ax,  ay,  az ,  which  satisfy  the  equation  (Morrill,  1961) 

a2  +  al  +  a2y  =  1.  (B.8) 

x  y  z 

Next  we  define  two  direction  cosine  matrices  [a,  ]  and  [  a, )  by  the  equation 
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Recall  that  the  4  by  4  matrix  [ y]  appearing  in  equation  (B.l)  is  defined  in  terms  of  the 
conductivity  tensor  [a]  by  equation  (20),  which  has  been  rewritten  in  equation  (B.10) 
below  for  convenience. 

(Yl  =  —  [o] .  (B.10) 

c 

The  complex  dielectric  tensor  [e]  (Born  and  Wolf,  1965)  is  defined  in  terms  of  the 
dielectric  tensor  [e]  and  the  conductivity  tensor  [a]  by  the  equation 

4k 

f el  =  [el  +  |o] .  (B.ll) 

Cl) 
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It  can  be  easily  shown  that  substitution  of  the  electromagnetic  plane-wave  solution  (B.2) 
back  into  the  Maxwell  field  equations  (B.l)  gives  the  following  matrix  equation 


cqe  cc2n 

/oi 

ICO 
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/oi 

a2  a,p 

f(12 

kc 
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/o2 

By  performing  the  matrix  multiplication  operations  suggested  in  equation  (B.l 2),  it  can 
be  easily  verified  that 

[aj  [/oil  =  [o]  and  [at]  [/02]  =  [oj  (B.13) 


and 

ICO  ICO 

[°4l  [/oil  =  [/02]  and  [a2]  [fm]  =  [/„,  1  •  (B.14) 

Equations  (B.13)  and  (B.14)  can  be  combined  to  form 

ICO  ICO 

[a,  +a2]  [/01]  =  \iyc  [/02]  and  [a,  +  a2]  [/„,]  =  e—  [/0]1 ,  (B.15) 

which  has  numerical  stability  attributes  not  found  in  equation  (B.14)  under  the  operation 
of  matrix  inversion.  It  can  be  easily  shown  that  the  sum  of  the  direction  cosine  matrices 
has  the  property 

[cq  +  cq]  [a,  +  a2l  =  -[/l  (B.16) 

where  [l]  is  the  4  by  4  identity  matrix.  This  implies  that  |oq  +  <xJ  has  a  multiplicative 
inverse  given  by 

[cq  +  a2] -1  =  -  [a,  +  a2] .  (B.17) 


Multiplying  either  equation  in  (B.15)  by  the  matrix  [oq  +  oqi  and  using  the  remaining 
equation  in  (B.15)  with  me  multiplicative  property  (B.16)  gives 


arpe 


k2c2 


=  1. 


(B.l  8) 


From  this  equation  we  obtain  the  following  expression  for  the  wavenumber  k 

c 


(B.19) 
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which  we  see  is  complex  for  a  conductive  medium.  With  the  use  of  equation  (B.19), 
equation  (B.15)  can  be  rewritten  in  the  finalized  symmetrical  form 

1/oiJ  =  ~‘Jj  lai  +  «2J  I/02]  and  I/02J  =  ~‘\j^  (a,  +  cx2l  I/01 J  -  (B-20> 

For  completeness,  equation  (B.20)  is  equivalent  to  the  following  vector  equations 


a  •Eg  =  0  and 

a.H0 

=  0 

(B.21) 

E„  =  -J[(axH„j 

and 

H>  = + 

(B.22) 

where 

E0=  (E„,  Eoy,  Eoz) 

and 

Ho~  (Hox’HorHoz). 

(B.23) 

Equations  (B.21)  and  (B.22),  and  hence  equation  (B.20),  express  the  transversal  nature  of 
the  electric  and  magnetic  fields  associated  with  the  monochromatic  electromagnetic 
plane-wave  solutions  (B.2). 

It  is  of  interest  to  note  that  the  complex  index  of  refraction  (Born  and  Wolf,  1965)  is 
given  by  the  equation 


In  addition,  the  complex  index  of  refraction  can  also  be  written  in  the  following  form 
(Born  and  Wolf,  1965) 

n  =  n  (1  +  ik)  (B.25) 


where  both  the  quantities  n  and  k  are  real.  The  quantity  n  is  the  ordinary  refractive 
index  and  k  is  the  attenuation  index.  The  quantities  n  and  k  may  be  expressed  in  terms 
of  the  material  constants  e,  |i,  and  c.  Setting  equation  (B.24)  equal  to  equation  (B.25) 
gives  the  relation 

,  ,  2n pa 

n2(1-K‘)  =  pe  and  n2K  =  .  (B.26) 


It  is  noted  that  experimental  values  of  the  optical  constants  n  and  k  for  a  variety  of 
materials  may  be  found,  for  example,  in  Born  and  Wolf  (1965). 
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APPENDIX  C 


ELECTROMAGNETIC  BOUNDARY  CONDITIONS 


The  electromagnetic  boundary  conditions  that  must  be  satisfied  at  the  planar 
interface  z  =  zp  separating  subregions  p  and  p+ 1  (see  figure  Cl  below)  are  contained  in 
the  matrix  equation  (33).  This  equation  has  been  rewritten  in  equation  (C.l)  below  for 
convenience 
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Figure  Cl.  Planar  interface  z  =  z  separating  subregions  p  and  p+1. 


The  matrices  fC,j  and  [C2]  appearing  in  equation  (C.l)  are  given  by  equation  (30). 
Equation  (C.l)  can  be  written  in  the  following  simplified  form 
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by  using  the  following  identity 
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The  matrices  [G,]  and  [G2]  are  given  by 
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(C.4) 


It  can  be  easily  shown  that  the  single  matrix  equation  (C.2)  is  equivalent  to  the  following 
four  boundary  conditions  in  vector  form  evaluated  at  z  -  zp: 

(1)  The  tangential  component  of  the  electric  field  is  continuous 

{zxE}p  =  {zxE}p  +  ,  (C.5) 


(2)  The  normal  component  of  the  electric  displacement  is  continuous 

{z«D}p  =  {z*D}p  +  ]  (C.6) 

(3)  The  tangential  component  of  the  magnetic  field  is  continuous 

{zxH}p  =  (:xH}p  +  ]  (C.7) 

(4)  The  normal  component  of  the  magnetic  induction  is  continuous 

{z-B}p  =  {z*B}p  +  1  (C.8) 


The  vectors  £,  H.D.and  B  appearing  in  equations  (C.5)  through  (C.8)  are  given  by 


£=  (E* 

E  E  ) 

and 

H  =  (Hx,HrHz) 

(C.9) 

D  =  e  ( Ex, 

E  ,  E  ) 

v’  z' 

and 

B  =  q  (Hx,  Hv.  Hz) . 

(C.10) 

The  unit  vector  z  points  along  the  positive  z-direction. 
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We  will  now  restrict  our  attention  to  equation  (C.2)  for  the  special  case  of 
monochromatic  plane-wave  solutions  previously  described  in  detail  in  appendix  B.  It  is 
a  well  known  fact  that  a  dielectric  multilayer  optical  medium  illuminated  by  a 
monochromatic  electromagnetic  plane-wave  will  give  rise  to  two  monochromatic  plane- 
waves  propagating  in  opposite  directions  within  each  layer  ot  the  multilayer  medium 
(Fowles,  1968).  The  electromagnetic  field  vector  associated  with  any  subregion  within 
the  multilayer  medium  may,  therefore,  be  represented  by  the  superposition  of  two 
electromagnetic  field  vectors 


X 

\ 

X 

i 

+ 

X 

h 

h 

R 

h 

(C.ll) 


Each  of  the  field  vectors  on  the  right-hand  side  of  equation  (C.ll)  has  the  monochromatic 
plane-wave  representation  given  by  equation  (B.2)  in  appendix  B,  that  is, 


X 

/oi 

h 

R 

fo2 

X 

foi 

h 

L 

/o2 

exp{i[r]r[k]R} 


exp{i[r\T[k]L} 


(CM) 


(C.13) 


The  subscripts  R  and  L  denote  wave  propagation  to  the  right  and  left  respectively. 
Substitution  of  equation  (C.ll)  into  equation  (C.2)  leads  to 


0 

I 

0 

G2 

j 

G,  0 

O  G, 


p+ 1 


+ 


p  + 1 


(C.14) 


For  equation  (C.14)  to  be  satisfied  for  all  values  of  time,  /,  the  angular  frequencies  of  the 
individual  plane-waves  in  the  two  subregions,  p  and  p+1,  must  satisfy  the  equation 

=  W,  =  {(0k}p+]  =  =  (C.15) 

Similarly,  for  equation  (C.14)  to  be  satisfied  for  all  values  of  *  and  y  on  the  planar 
boundary  z  =  zp ,  the  at-  and  y-components  of  the  wave  vectors  associated  with  each  of 
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the  individual  plane-waves  in  the  two  subregions,  p  and  J,  must  satisfy  the  equations 


0,0,= 

0,0,= 

0,0, „ 

=  O,0r„  =  K 

(C.16) 

0,0,= 

0,0,., 

=  0,0 =  V 

(C.17) 

From  equation  (B.19)  it  is  apparent  that  the  wavenumbers  associated  with  the  R  and  L 
waves  in  any  subregion  within  the  multilayer  medium  must  be  equal,  namely 

{**>,  =  (*l>,  =  kp  and  =  U'L}p+1  =  (C.18) 

With  the  use  of  equations  (C.16),  (C.17),  (C.18)  and  (B.7),  (B.8)  from  appendix  B,  it  can  be 
easily  shown  that  the  z-components  of  the  wave  vectors  associated  with  each  of  the 
individual  plane-waves  in  the  two  subregions,  p  and  p+1,  must  satisfy  the  equations 


(Kr}  p  ~ 

-o,op 

-  00, 

(C.19) 

{kzR}p  +  ]  = 

-0,0,., 

=  00 

(C.20) 

With  the  use  of  equations  (C.12)  through  (C.20),  the  matrix  equation  (04)  takes  on  the 
finalized  form 
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exp  ( -izpkz )  = 


(C.21) 
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p+  1 
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